Pseudogap in high-temperature superconductors from realistic Frohlich and Coulomb 

interactions 
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It has been recently shown that the competition between unscreened Coulomb and Frohlich 
electron-phonon interactions can be described in terms of a short-range spin exchange J v and an 
effective on-site interaction U in the framework of the polaronic t-J v -U model. This model, that 
provides an explanation for high temperature superconductivity in terms of Bose-Einstein conden- 
sation (BEC) of small and light bipolarons, is now studied as a charged Bose-Fermi mixture. Within 
this approximation, we show that a gap between bipolaron and unpaired polaron bands results in a 
strong suppression of low-temperature spin susceptibility, specific heat and tunneling conductance, 
signaling the presence of normal state pseudogap without any assumptions on preexisting orders or 
broken symmetries in the normal state of the model. 

PACS numbers: 74.20.-z, 74.72.-h, 71.38.-k 



Experimental evidence of finite charge/spin pseudo- 
gap (PG), namely a depression of the electronic ex- 
citation spectrum at a temperature T* well above 
the critical temperature T c . have been widely advo- 
cated as one of the most significant signatures of hid- 
den orders or broken symmetries in the underdoped 
regime of high-temperature superconductors. Although 
a plethora of different techniques have been used for 
the investigation of the PG, such as tunneling mea- 
surements, angle-resolved photoemission spectroscopy 
(ARPES), nuclear magnetic relaxation (NMR), Raman 
and neutron scattering '", a microscopic theoretical de- 
scription is still missing . 

Early NMR measurements, in which spin excitations 
are probed, revealed an anomalous depression in the tem- 
perature dependence of the Knight shift for underdoped 
samples of YBCO J . Such a depression, well-known in 
BCS-like superconductors as induced by the formation 
of a spin-singlet state, allows some researcher to main- 
tain that PG features could be explained in terms of pre- 
formed Cooper pairs ' 5 or as signatures of the suppression 
of low-energy antiferromagnetic fluctuations How- 
ever, a later interpretation of NMR data as a consequence 
of a depression in the electron density of states (DOS), 
results in a number of controversial debates mainly based 
on the awareness that, in principle, any instability (e.g. 
stripes, charge/spin density waves, polaron formation) 
might result in an energy gap. Therefore, the conclu- 
sion that pseudogap features might not necessarily im- 
ply spin-singlet formation, paves the way to other theo- 
ries and phenomcnological models in which PG emerges 
as a consequence of SU(2) rotation 10 , coexistence of 
charge and spin density waves 11 , inhomogeneous charge 
distributions 12 ' 13 . 

Recently, strong evidence of particle-hole symmetry 
breaking in the pseudogap state of Bi2201 shed doubts 
about the possibility to consider the pseudogap as a 



precursor of a Cooper pairing superconducting gap in 
the normal state 14 . Also some earlier and more recent 
ARPES experiments 15 emphasized the role of electron- 
lattice coupling as an unavoidable ingredient for the char- 
acterization of both normal and superconducting state 
of high-tempcrature superconductors. Actually, on the 
phenomcnological level, the pseudogap was originally ex- 
plained as half of the bipolaron binding energy . 

In this context, we report our study on the pseudogap 
in the framework of the microscopic t-J p -U model: 

H = tij6<r<r'C\Cj + "-mt«m| + 

i,j m 

+2 J2 MU, m - n) ( S m • S n + -n m n n J (1) 

that accounts for realistic Coulomb repulsion and strong 
electron-phonon (Frohlich) interaction in terms of an ex- 
change coupling J p ' and a residual on-site correla- 
tion U which, limiting the double occupancy, reduces 
J p from its bare value 19 . Here Cj, c\ are polaron an- 
nihilation and creation operators where i = (m, a) and 
j = (n, a 1 ) include both site (m, n) and spin (cr,cr') in- 
dices and the sum over n/m counts each pair once only. 
n m = "mt + "mi, and n mt ^ = c mt ^c m1 -4 are site oc- 
cupation operators and S m = (1/2) J2<t.<j> c lna~^ °<y' c ma> 
is the spin 1/2 operator (r are the Pauli matrices). It 
is worth recalling that our t-J p -U model describes carri- 
ers doped into the charge-transfer Mott-Hubbard (or any 
polar) insulator, rather than the insulator itself, different 
from the conventional Hubbard U or t-J models. The 
bare Hubbard-C/ on the oxygen orbitals (where doped 
holes reside) in a rigid cuprate lattice is of the same or- 
der of magnitude as the on-site attraction induced by 
the Frohlich EPI (ss leV to 2cV) so that the residual 
Hubbard U could be as large as a few hundred meV 19 . 



2 



Hereafter we restrict our analysis to nearest-neighbor J 
and t on a square lattice. 

We show that, without any ad- hoc assumption on the 
relative strength and the range of Coulomb and electron- 
phonon interactions, PG features naturally appear as a 
consequence of a thermal-induced mixture of polarons 
and bipolarons, providing a non-Fermi liquid description 
for the normal state. 

Early numerical and analytical studies on the t-J p and 
t-Jp-U models pointed out that the ground state of the 
system can be described as a coherent supcrfluid of inter- 
site small bipolarons 18,19 . As follows from a straightfor- 
ward calculation of the static (t = 0) ground state config- 
uration, bipolarons repel each other due to the presence 
of a short-range effective bipolaron-bipolaron repulsion 
Ef,b = (2 — y/3)J p (U). Furthermore, they do not attract 
single polarons. Hence there is no tendency to cluster- 
ing. At this point, it is convenient to apply a potential 
shift to the Hamiltonian: H — >• % — \Eq J2m n m- Here 
Eq is the two-particle ground state energy; for small t, 
Eq = —J p. Because of the aforementioned bipolaron- 
bipolaron repulsion, all energy eigenstates of the shifted 
Hamiltonian are non-negative. The shift has no physical 
effect, as it will be absorbed into the chemical potential, 
but provides a more intuitive visualization of the spec- 
trum (Fig.l). 

As long as we are in the low-density regime, bipo- 
laron and unpaired polaron interactions are negligible 
and the model can be described in terms of an ideal (non- 
interacting) Bose-Fermi mixture of bipolarons (bosons) 
and unpaired polarons (fermions) in which bipolaron and 
unpaired polaron densities n p ^{T) are related to the dop- 
ing x <C 1 as 20 : 

x = 2n 6 (T, fi b ) + n p {T, fi p ) . (2) 

Here ji p> b represents the polaronic/bipolaronic chemical 
potential. Charge conservation allows us to fix the chem- 
ical potential \x of the whole mixture. In particular, ac- 
cording to detailed equilibrium and total energy conser- 
vation one readily obtains \i = \i p = /Uft/2. Then the 
thermodynamic properties of the Fermi-Bose mixture at 
equilibrium can be easily derived in terms of the ther- 
modynamic potential f2(a;,T), defined as a function of 
temperature T and chemical potential /i 21 . 

In the dilute limit all the thermodynamic quantities 
enjoy the additivity property therefore £l(x,T) can be 
expressed as the sum of bipolaron/unpaired polarons po- 
tentials £lb,p{nb,p,T) with: 



fi P ,b = TfeT / deJ\f P}b (e) In li exp 



fJ-p,b 



k B T 



(3) 

Here A4, P (e) is the density of states in the bipo- 
laron/unpaired polaron band, respectively. For the sake 
of simplicity, hereafter the DOS of the model is approxi- 
mated as: 



\Afp(e) = A p Q (e - A) 9 (A + 2w p - e) 
\M b (e) = A b Q(e)Q(2w b - e) 



(4) 



In a square lattice with nearest-neighbour J pi there are 
two bound bosonic states per site, corresponding to the 
two bonds per site and two fermionic states per site due 
to spin degeneracy. Thus rectangular bands (including 
the degeneracy factors) need to satisfy the constraint 
J J\f Pj b(e)de = 2 from which we have A P: b = l/w Pi b- Here 
Ap^b and w Pj b are intensity and half-bandwidth of po- 
laron and bipolaron terms, respectively, A represents the 
binding energy per polaron. 
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Figure 1: (Color online) Left panel: bipolaron (bottom) and 
unpaired polaron (top) bands with the corresponding DOS: 
N p , b (E) (filled area) for U / J P {U) = 2.0 and t/J p {U) = 0.1. 
Here Eq = 2Eq — Eq, where Eq is the two-particle ground 
state energy. Dashed lines on the DOS represent the resulting 
Heaviside-theta approximation. Right panels: DOS parame- 
ters versus t/J p (U) for different values of the ratio U /J P (U). 

Importantly, DOS and model parameters are linked 
together: A w J p (U)/2 at t/J p <C 1 while bandwidths 
Wp.b are related to the ratio t/J p (U) as reported in Fig.l. 
Such an approximation is able to get an insight in a 
qualitative microscopic description of the pscudogap. It 
allows to obtain analytical expressions for the relevant 
physical properties. In particular, bipolaron and polaron 
densities, calculated as dd\f p ,b{c)f P ,b(c, T), where 
f p ,b(e,T) = [exp (e - fi Ptb )/k B T ± 1] _1 is the Fermi- 
Dirac/Bose-Einstein distribution function, are expressed 
as: 



n b (T) 



n p (T) 



1 + *2lm 



sinh 



sinh ( — 



k B T 
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cosh(^ii) 



cosh 



\2k B T ) _ 



(5) 



from which /i can be calculated sclf-consistcntly accord- 
ing to Eq.2. 

Our results on chemical potential and particle density, 
reported in Fig. 2 for a fixed total number of particle x, 
show that different temperature behaviors arise in po- 
laron and bipolaron densities depending on the value 
of A and on the competition between the pairing in- 
teraction J P (U) and hopping term (Fig.l). In particu- 
lar, as follows from the left panel of Fig. 2, the bipolaron 
density decreases with increasing temperature resulting 
in a crossover at T = T* when half of the bipolarons 
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Figure 2: (Color online) Left panel: relative bipolaron 
(solid) /unpaired polaron (dashed) density versus temperature 
for different values of the gap A. The dotted line represents 
the total particle density x = 2n b + n p . Right panel: linear 
dependence of the ratio A/T* with respect to ln(l/x) for dif- 
ferent value of the gap A. In the inset the doping dependence 
of T*/A (symbols) is compared with the exact analytical de- 
pendence (line) obtained in the zero-bandwidth limit. Here 
T* is the crossover temperature at which nt(T*) = 2n p (T*). 



are dissociated and the charge is equally distributed be- 
tween polarons and bipolaron (rib(T*) = 2n p (T*)). As 
shown in the right panel of Fig. 2 the ratio A/fc^T* varies 
linearly with ]n(l/x) in a wide range of doping, with 
k B T* = 2A/ln((4/a- - 1) 2 /(1 + 8/x)) for w p , b /A -> 
in agreement with exact analytical calculations in the 
narrow-band limit. 

The thermal-induced population of the unpaired po- 
laron band results in a number of anomalous features 
that can be observed in specific heat C(T) and static uni- 
form spin susceptibility x s (T,h) (h is an external mag- 
netic field). Recalling that C(T) = d(E)/dT, where 
(E) = (E p ) + (Eb) is the total energy of the system 
with (E Ptb ) = JdeEAf p , b (e)f p ,b(e,T), and X s(T,h) = 
—d 2 £l/dh 2 , we have: 



Xs(h,T) = 



2 

2k b T 



dej\f p (e) 



C{T) = k B T A k 

k — b,p 

Here fig is the Bohr magneton and 
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while I pb \ x ) = / dx ^±\yi 1S expressed in terms of the 
polylogarithm function Li s (z) — X^feLi z fe /fc s as: 



r (»), v _ I li^Tlnaie 1 ) 

1 p,b \ X ) ~ \ „ ( le 



T 2 In (I ± e 1 ) ] =p Li2 fre 



n = 1 
n = 2 



. (9) 

As follows from Fig.l, for t/J p (U) w 1 and for J7 > 1 
{J P (U — > oo) — > 19 ), the gap A goes to zero. In this case 
bipolaron and unpaired polaron bands are completely 
overlapped with a single peak in the specific heat coeffi- 
cient 7 (T) = C(T) jT leading to -y{T) cx 1 /T, Fig.3. Con- 
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Figure 3: (Color online) Specific heat coefficient y(T) = 
C(T) /T versus temperature plotted for different values of gap 
A (left panel) and doping x (right panel). Here Wb. p is the 
half-bandwidth of the bipolaron/unpaired polaron band. 



sistently, the paramagnetic response of the system repro- 
duces the standard Curie law with Xs(T, h = 0) oc 1/T 
(Fig. 4) since fermions are non-degenerate under the con- 
ditions here. 

On the contrary, in the opposite regime the presence 
of a finite gap results in a non monotonic j(T) depen- 
dence induced by the superposition of two main peaks 
due to intra-band and bipolaron to unpaired polarons ex- 
citations. While any finite temperature can induce intra- 
band excitations, bipolaron dissociation requires temper- 
atures of the order of the gap. Therefore the separation 
between the two peaks increases with increasing A and 
results in a strong suppression of the 7 (T) coefficient in 
the region in which the intensity of the intra-band peak 
falls off (ksT sa Wb). Specific heat and static uniform 
spin susceptibility are both determined by the average 
density of electronic states therefore, as one would ex- 
pect, the same tendency is also observed in the para- 
magnetic response function Xs(T,h = 0). In fact, as 
confirmed in Fig. 4, x(T, h = 0) drops to zero in the low 
temperature regime in which the population of the bipo- 
laronic (singlet) band has its maximum and follows the 
standard Curie law, already described in the A = case, 
in the high-temperature regime in which the population 
of the polaronic band become dominant. The same fea- 
tures also appears in the presence of an external mag- 
netic field until /xg/i ~ A. For hb h > A the magnetic 
field induces a finite magnetization in the system with a 
singlet/triplet phase transition at fish = A signalled by 
a discontinuity in the spin susceptibility at T = 0. 

Finally, let us discuss the tunneling conductance of a 
normal mctal-bipolaronic superconductor (NS) junction: 



ctns(V) 



dI NS {V) 
dV 



(10) 



in which V is the bias and Ins{V) is the current 
flowing through the junction. According to the the- 
ory of extrinsic and intrinsic tunneling in bosonic 
superconductors 24,25 , Ins(V) can be calculated starting 
from the following tunneling Hamiltonian: 



H„ e =P 



(11) 
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Figure 4: (Color online) Spin susceptibility (Eq.6) versus tem- 
perature plotted for different values of the external magnetic 
field h. Here A is the gap between bipolaron and unpaired po- 
laron bands, w p is the half-bandwidth of the polaronic band. 



that accounts for single particle tunneling processes. 
Here c v , p\, and 6^ describe the annihilation of a car- 
rier in the metallic tip in state v and the creation of 
a single polaron or a composed boson in the supercon- 
ductor in state v' or r\ respectively, N is the number 
of lattice cells. P and B are tunneling matrix elements 
respectively with and without the involvement of a bipo- 
laron (generally B > P ). The tunneling current is 
Ins(V) = eiWn^s — Ws^n) where Wx^y represents 
the tunneling probability of transition, per unit time, 
from the X to the Y side of the junction. According 
to the Fermi golden rule we have: 



)=^A m A p l [ \P\ 2 i 



[M£'-eV0-MO] 



r2w b rA + 2w p 

+ \B\ 2 A b / d V / d?f F (v - C - eV)f F {0 

JO J A 

r2w b rA + 2w p . 



(12) 



where /f(£) = [cxp (^/fcgT) + 1] 1 is the Fermi distribu- 
tion functions associated to normal metal and polaronic 
band; /s(?y) = [cxp (i]/^T) — 1] _ represents the Bose 
distribution function associated to the bipolaronic band. 
Here we have used a constant DOS for the normal metal 
with J\f m (e) = A m . 

As reported in Fig. 5, in the same regime in which pseu- 
dogap features arise in the specific heat (Fig. 3) and spin 
susceptibility (Fig. 4), our results on the tunneling con- 
ductivity confirm a strong depression of u{V) at zero bias 
for ksT < A. Importantly, despite the lack of van Hove 
singularities in the DOS (Fig.l), our data perfectly repro- 
duce the asymmetry between negative and positive bias 
conductance providing a further confirmation that the 
van Hove singularity, not observed in many experiments 
such as momentum integrated photocmission- 1 ', does not 
play any role in the tunneling. Our model allows us to 
describe the doping dependence of the asymmetry coef- 
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Figure 5: (Color online) Left panel: density plot of the 
normalized conductivity <jnm{V) /onm{A) in the ksT/A - 
(eV — fi(T))/A plane. Right panel: doping dependence of 
the asymmetry coefficient R(x,T), Eq.f3. Numerical re- 
sults obtained by integrating the normalized conductivity 
a(eV) /a(A) from to ±A for different values of the tem- 
perature are compared with experimental results in cuprates 
(from Ref.25). 



ficient R(x, T) defined as: 

S°_ A a(eV)dV _ I NS (-A) 



R{x,T) 



J*a{eV)dV Ins(A) 



(13) 



We recall that with increasing doping the bipolaron den- 
sity increases while the unpaired polaron density remains 
almost zero for ksT < Wb (see Fig. 2, left panel). In this 
regime the polaronic contribution to the tunneling con- 
ductance remains constant while the bipolaron contribu- 
tion, different from zero only in the positive bias regime, 
scales linearly with the bipolaron density as clearly fol- 
lows from Eq.12 if one neglects the bipolaron energy 
dispersion in the narrow bipolaron-band limit. As re- 
ported in Fig. 5, this explanation is also supported by 
a good agreement between our data and experimental 
measurements of the asymmetry coefficient R(x, T) in a 
wide range of cuprates superconductors. Importantly, 
it is worth noting that numerical data for R(x, T) have 
been calculated by integrating the normalized conductiv- 
ity <r(eV) /<t(A), therefore do not depend on the partic- 
ular choice of the tunneling matrix elements B, and P 
in Eq.12. As one would expect, the only relevant quanti- 
ties are t, J p {U), U that, unlike in other theories, in the 
t-Jp-U model can be fixed from the measurable material 
properties 17-19 . 

In conclusion, we have described the normal state of 
the polaronic t-J p -U model as an ideal Bose-Fermi mix- 
ture in the low density limit. By approximating the DOS 
of the model as rectangular functions (4), we have pro- 
vided analytical expressions for bipolaron and unpaired 
polaron densities (5) and most of the relevant response 
functions such as specific heat (7), spin susceptibility 
(6). Our analysis pointed out that in the presence of 
a finite gap A between bipolaron and unpaired polaron 
bands, the model exhibits remarkable features of pseudo- 
gap opening signaled by a depression of specific heat, spin 
susceptibility and tunneling conductance. Furthermore, 
as the result of the screening 22 ' 2 ' 5 , the pseudogap A falls 
with doping through the dielectric function so that the 
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crossover temperature of our model falls as well. Impor- 
tantly, different from any other theories proposed so far, 
pscudogap features naturally appear as a consequence of 
a thermal-induced mixture of polarons and bipolarons, 
without any ad-hoc assumption on relative strength and 
the range of Coulomb and electron-phonon interactions 



or preexisting orders. 
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